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Abstract 

The theory of phase ordering kinetics for the 0(2) model using the gaussian 

auxiliary field approach is reexamined from two points of view. The effects of 

fluctuations about the ordering field are included and we organize the theory 

such that the auxiliary field correlation function is analytic in the short-scaled 

distance (x) expansion. These two points are connected and we find in the 

refined theory that the divergence at the origin in the defect-defect correlation 

function g(x) obtained in the original theory is removed. Modifications to the 

order-parameter autocorrelation exponent A are computed. 
PACS numbers: 05.70.Ln, 64.60. Cn, 64.60.My, 64.75.+g 
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I. INTRODUCTION 



The phase ordering kinetics of systems with continuous symmetry, such as the 0(n) 
model, is particularly interesting because of the topological defect structures produced: vor- 
tices and strings for n = 2 and monopoles for n = 3 |. While much is understood about 
the theory of growth kinetics for the 0(n) model, there are some interesting unresolved 
problems associated with the short-distance behaviour of the defect-defect correlation func- 
tion g(x). In the theory developed in one expresses the order parameter ip(R,t) as a 
local non-linear function of an auxiliary field m(R, t) which is physically interpreted as the 
distance, at time t, from position R to the closest defect. One of the physical motivations 
for introducing m(R, t) is that it is smoother than the order-parameter field. Sharp inter- 
faces or well defined defects produce a non-analytic structure in the order-parameter scaling 
function T{x) at small-scaled distances x which, physically, is responsible for the Porod's 
law decay seen scattering experiments 0. The expectation, however, is that the auxiliary 
field correlation function f(x) will be analytic in this same distance range. In the case of 
a scalar order-parameter these expectations are supported by the theory ||. However for 
n > 1, as pointed out in 0, this is not the case. One finds a weak non- analytic component 
in / and, more significantly, for n = 2 one can trace this non-analytic component to an 
unphysical divergence in g(x) at small x. |J. This divergence is not seen in simulations |3j 
or experiments @ for n = d = 2 where g(x) apparently approaches zero at the origin. 

In this paper we focus on the case n = 2 and show how these problems can be resolved 
by taking seriously the assumption that the correlations of the auxiliary field are indeed 
smoother than those of the order parameter. We find that it is possible to rearrange the 
theory such that / is analytic in x if we extend the theory to include fluctuations about the 
ordering field and treat the separation between the ordering field and the fluctuation fields 
carefully. This is accomplished by introducing a new field O which is constructed to ensure 
that the fluctuations are small, while at the same time compensating for the non-analyticities 
in /. It is important to note that we work at zero-temperature, so the fluctuations are not 
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thermally driven ||. Rather, we will see that the correlations in the fluctuations are slaved 
to the correlations in the order parameter. The theory of Ohta, Jasnow and Kawasaki (OJK) 
H basically avoids this entire discussion by simply assuming that / has a gaussian form |jTD[ . 
There is no self-consistent determination of / in that theory. 

It is well established that for late times following a quench from the disordered to the 
ordered phase the dynamics obey scaling and the system can be described in terms of a 
single growing length L(t), which is characteristic of the spacing between defects. In this 
scaling regime the order-parameter correlation function 

(7(12) = <#1) • <?(2)) (1.1) 

has a universal equal-time scaling form 

C{12) = $F(x), (1.2) 

where ipo is the magnitude i[> = of the order-parameter in the ordered phase. Here 
we use the short-hand notation where 1 denotes (Ri,ti), and define the scaled length x as 
x = R/L(t) with R = |R| = |R 2 — Ri|. It is also well established that, in the scaling regime, 
L(t) ~ P* where t is the time after the quench. For the non-conserved models considered 
here the exponent <ft = 1/2. Another measurable quantity is the exponent A governing the 
decay of order-parameter autocorrelations, and defined by 

C(0 ) i,f)~-i-fori»f, (1.3) 

This non-trivial exponent can be computed theoretically, along with the scaling function T 



Q. The predictions for A are in excellent agreement with simulation results Jll|JT2| . The 
theoretical predictions for T are also in good agreement with simulations @ . 

The dynamics of the defect structures themselves is amenable to theoretical treatment 
||. In this paper we shall mainly be interested in the case n = d, where the defects are 
points. For n — d — 1 the defects are strings, but the analysis follows closely that for point 
defects and yields qualitatively similar results. The density of point defects, for n = d, is 
defined as 



p(l) = J £q a S{R 1 -x a {t 1 )) (1.4) 

a 

where x Q ,(t 1 ) is the position at time t\ of the ath point defect, which has a charge q a . 
Correlations in p, 

G(12) = (p(l)p(2)>, (1.5) 

at equal-times ti = t 2 = t can be shown || to decompose into two parts 

G(R,t) = n Q (t)5(R) + g(R,t). (1.6) 

The first term n (t) represents defect self-correlations and is just the total unsigned number 
density of defects at time t. We will be primarily concerned here with the second term 
g(R,t) which measures the correlations between different defects. In the scaling regime it 
can be shown || that n (t) ~ L~ n {t) and that g(R,t) has the form 

g(K,t) = J ^ r) g(x). (1.7) 

where g(x) is a universal scaling function. 

In the next section we present the 0(n) model and describe the mapping between the 
auxiliary field and the order parameter. In Section III we discuss the separation of the 
equation of motion into an equation for the evolution of the ordering field and an equation 
for the dynamics of the fluctuations. The main analytical results of the paper are presented 
in Section IV, where we discuss how the quantities J 7 , A and g are determined through 
the solution of a non-linear eigenvalue problem. In Section V we calculate the correlations 
in the fluctuations, assuming that the fluctuation field u and the auxiliary field m form a 
set of coupled gaussian variables. Our numerical analysis of the new non-linear eigenvalue 
problem is presented in Section VI and the results are discussed in the final section, which 
also addresses more general issues that indicate directions for future research. 



II. MODEL 



We consider the 0(n) model, which describes the dynamics of a non-conserved, n- 
component order-parameter field ^(1) = (V'i(l), - • • , ip n (l)). To begin we will work with 
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general n; however, later we will focus on the interesting case n = 2. As in previous work in 
this area ||], the dynamics are modeled using a time-dependent Ginzburg-Landau equation 

dt ~ $ ■ [ZA) 

We assume that the quench is to zero temperature where the usual noise term on the right- 
hand side is zero []13| |. T is a kinetic coefficient and F[ip] is the free-energy, assumed to be 



of the form 



FM= J d d r(^\ 2 + Vty)) (2.2) 

where the potential V[ip] is chosen to have 0(n) symmetry and a degenerate ground state 
with if> = ip . Since only these properties of V will be important in what follows we need 



not be more specific in our choice for V JTJ]. With a suitable redefinition of the time and 



space scales the coefficients T and c can be set to one and (|2.1|) can be written as 

i = vv '-^f- (2 - 3) 

It is believed that our final results are independent of the exact nature of the initial state, 
provided it is a disordered state. 

The evolution induced by (|2.3|) causes ip to order and assume a distribution that is far 
from gaussian. It is by now standard to introduce a mapping between the physical field ij) 
and an auxiliary field m with more tractable statistics. We can decompose ip exactly as 

$ = a[m] + u. (2.4) 

The utility of this decomposition lies in our ability to create a consistent theory with the 
mapping a chosen to reflect the defect structure in the problem, and u constructed to be 
small at late-times. Thus u represents fluctuations about the ordering field a. The precise 
statistics of the fields m and u will be specified below. 

The defect structure [|T5[] is naturally incorporated by using the Euler-Lagrange equation 
for the order-parameter around a static defect in equilibrium, 



W\ = ^1, (2.5) 

to determine the functional dependence of a on fh. The defects are then the non-uniform 
solutions of ( p.5| ) which match on to the uniform solution at infinity. Since we expect only the 
lowest-energy defects, having unit topological charge, will survive to late-times the relevant 
solutions to ( p.5| ) will be of the form 

<j[fh] = A(m)fh (2-6) 

where m = \fh\ and fh = fh/m. Thus the interpretation of fh is that its magnitude represents 
the distance away from a defect core and its orientation indicates the direction to the defect 
core. We expect m, away from the defect cores, to grow as L in the late-time scaling regime. 
Inserting ( |2.6|) into (|2.5| ) gives an equation for A, 

71 — 1 

VIA-— -A-V'[A]=0, (2.7) 

fit 

where the prime indicates a derivative with respect to A. The boundary conditions are 
A(0) = 0, A(oo) = if) . An analysis of ( |2.7| ) for n > 1 and large m yields 



A(m) = ip 



K 

+ ■ ■ 



(2.8) 



where k — (n—1) /V"[ipo] > 0. The algebraic relaxation of the order-parameter to its ordered 
value is a distinct feature of the 0(n) model for n > 1. In the scalar case (n = I) ip relaxes 
exponentially to ipo away from the defects. 

We shall also be interested in the stability matrix defined by 

™-£g <-> 

VIA] 

= v"[A]& t a :i + —^(5 iJ -a i a J ). 

By the definition of ipo and what we mean by equilibrium we have 

V'[i[> ] = (2.10) 
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and 

l/'Vo] = ql > 0. (2.11) 

These results give 

WijbPo] = qfa&j (2.12) 

which is purely longitudinal. This reflects the fact that, in equilibrium, the longitudinal 
fluctuations have a "mass" q$ while the transverse fluctuations, or spin-waves, are massless. 

III. SEPARATION OF EQUATIONS OF MOTION 

In this section we develop the equations of motion satisfied by the fields a and u. Let us 
first define 

and rewrite the equation of motion (|2.3| ) in the form 

^(ff + u) = W\a + u)- H\a + u}. (3.2) 

We then quite generally assume that a satisfies the equation of motion 

^ = VV-V^ + (3.3) 

where is as yet unspecified. Clearly, subtracting this equation from the equation of motion 
( jOD and using (|2.5|) we obtain 

^ = V 2 u-H[a + u] +H[a] -6. (3.4) 

To this point things are quite general since we have not specified 0. A key point is that 
must be chosen such that u does indeed represent a fluctuation. This means that in the 
scaling regime we can treat u as small and keep only leading powers of u in the equations 
of motion for a and u. Equation ( |3.4|) , to leading order, is then given by 



(I'll ■ 

-^- = V 2 u i -W ij [a]u j -e i (3.5) 

where a sum over the index j is assumed. 

We now assume that is a function of m only. This means that a satisfies a closed 
equation, while u is slaved by m. We will choose the form for so that the correlation 
function / for in is analytic for short-scaled distances. As we shall see this is a rather 
constrained process. 



IV. ANALYSIS OF THE a DEGREES OF FREEDOM 

A. Construction of 6 

If we set G equal to zero in fl3.3| ) we obtain the equation used previously to determine 
the a correlations @]. This choice decouples a and u. The equation for u would then 
separate into a (massless |Tj|) diffusion equation for the transverse piece ut and an equation 
for the longitudinal piece ul with a mass term —q^UL- However, the equation for a would 
necessarily lead to non-analytic behaviour in / at short-scaled distances and ultimately to 
an unphysical divergence in g(x) at small x. We must choose G so that f(x) is analytic for 
small x. The form we can use for G is determined by the following observations: 

(i) Q must be odd under m — > — m. 

(ii) Q must scale as 0(L~ 2 ) in the scaling regime if it is to compensate for the terms in 
the equation of motion which lead to the non-analyticities in /. This will also allow us to 
treat fluctuation since it will imply u ~ L . 

It is not easy to construct a variety of functions of m which are independent and satisfy 
(i) and (ii). We propose the general form 

e^^+E^Vmf^v^ (4.i) 

where [Vmf = Yh=i E"=iPi m a] and all of the cue, i > are assumed to be of 0(1). One 
can think of including other quantities like (ip^ — a 2 )a but these, in the scaling regime, are 
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equivalent to V^er. It is interesting to note, using the definition ( |2.5| ) for V^cr, that is 
longitudinal. 

For the purposes of this paper we will only consider constructing f(x) to be analytic 
through terms of C(x 4 ). To satisfy this requirement it is sufficient to set £ max = 2 in Q4.1|) . 
The equation for a (|3.3| ) is then of the form 

B = (4.2) 

where we define, for later convenience, 

B = 8 t a- V 2 a + V 2 m a - -^-a - u{7 2 m ff - uu 2 [V m] 2 V 2 m a . (4.3) 



B. The Gaussian Approximation 

To complete the definition of the model one must specify the form of the probability 
distribution for the auxiliary field fh. Forcing a to satisfy the exact equation of motion ( |f ,2| ) 
is tantamount to solving the problem exactly, and will determine a probability distribution 
for m which is complicated and extremely difficult for purposes of computation. Progress 
can be made if one imposes the weaker constraint 

(B(l) ■ a{2)) = 0. (4.4) 

This equation allows one to insure that B(l) is reasonably small at late-times but gives one 
the flexibility to choose a suitable probability distribution. The simplest choice is a gaussian 
probability distribution for m with the correlation function Co(12) explicitly defined through 

(m i {l)m j (2)) = 6 ij C Q {12). (4.5) 

The system is assumed to be statistically isotropic and homogeneous so C*o(12) is invariant 
under interchange of its spatial indices. For future reference we also define the one-point 
correlation function 
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Sn(l) = Cn(ll) 



(4.6) 



and the normalized correlation function 



/(12) 



Co(12) 
5 (12) 



(4.7) 



with 5o(12) = ySo(l)So(2). As discussed above it is expected that both Co and So grow as 
L 2 at late times. The gaussian approximation, which has been successful in describing the 
correlations in these systems, forms the basis of almost all present analytical treatments of 
phase-ordering problems [0,H. Efforts to go beyond to gaussian approximation are defined 



in E7.1 



The functional dependence of C(12) and g(x) on / can be derived without reference 
to the dynamics contained in (|4.4j) . Using ( |2.4| ). (|2.6|) , and ( |2.8| ) C(12) can be written to 
leading order in l/L: 



C(12) = $(m(l)-m(2)) 
= ^(12). 



(4.8) 



Assuming gaussian statistics for fa we obtain [|4|,|T9f 



H12) = ^fi 2 



1 n + 1 



2' 2 



1 1 n + 2 



2'2' 2 ' 



/ 2 (12) 



(4.9) 



where B is the beta function and F is the hypergeometric function. Within the gaussian 
theory, g(x) is given by [|] 



g[x) = n\ 



n-l 



dh 
dx 



(4.10) 



with h = — 7/'/27r and 7 = 1/a/1 — p. The defect density is given by 



n (t) 



2 n W 2 r(l + n/2) 



S 



(2) 



nS (t) 



re/2 



(4.11) 



with 



^ 2) = -([Vmf). 



1 



(4.12) 
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C. Order Parameter Correlations 



With the specification of the probability distribution for m the constraint (O) leads to 



an equation that allows one to compute correlations in the order parameter a. The quantity 
S K o> ( 1~T2 ), which will later appear in the definition of the length scale and in the formula 



for the autocorrelation exponent A, is determined through condition (f4.4|), with 2 = 1: 



\d tl (a\l)) - {a{l) ■ V?*(l)> - j^&Q)) + (1 " <*X*(1) " V^(l)) 

-^ 2 (a(l)-[V 1 m(l)] 2 V^(l))=0. (4.13) 



Equation ( |4 . 1 3| ) can be simplified by using the following identities, true for gaussian averages, 

(a(l) ■ V?<r(l)> = -4 2) (V m • [5(1) • V m a(l)]> + Sf\a{l) • V^x(l)> (4.14) 
(<x(l) ■ [V 1 m(l)] 2 V^a(l)) = nSi 2) (a(l) • V 2 m <r(l)>, (4.15) 

and observing that, at late-times, the dominant term in ([4.13|) is 

<*(1) • V^(l)) = -^y In S (l) forn = 2 

^ n_1 forn>2. (4.16) 



S (l) n- 2 

Evaluating ( |4.13| ) to leading order in 1/L one has 



= 1 ; (4- 17) 

1 — 2fi(n — 2)u)o/n(n — 1) + ncu2 



where we have defined the scaling length 



7T 



5 (t) 



L\t) = -^ = At. (4.18) 

Note that for n = 2 the term with u does not appear in Sq because it is dominated by 
the 0{L~ 2 In L) terms in Q4.13j ). There is a further simplification of ([4.171) for n = 2 since 
later we will have to set iO\ + 2Sq lu 2 = to ensure that the correlations in u remain finite. 
With this relation we have 

= 1 for n = 2, (4.19) 
11 



(2) ~* 

which is the value for Sq obtained previously for all n when = 0. 

Equation ( |4.4|) directly determines the time evolution of the two-point order-parameter 
correlations and is given explicitly by: 

u 



L 2 (ti 



(a(l)-a(2)) +(1-^)^(1)^(2)) 

-^ 2 (V ? 2 n a(l)[V 1 m(l)] 2 -a(2))=0. 



(4.20) 



Equation (|OoD can be re-expressed as an equation for the order-parameter correlation 
function T (|4.9|) by means of the following identities: 

(V^(l)[V!ni(l)] 2 ■ 9(2)) = nS^(V 2 m 9(l) ■ 9(2)) + [V 1 C (12)] 2 (V 2 „a(l) ■ V 2 m a(2)) (4-21) 



(Vi?(l) " *(2)> 



> 2 



(V^cr(l) • V 2 n a(2)) = Z»(f df F + fdjF) 

On 



(4.22) 
(4.23) 



where we use the shorthand notation / = /(12), T = J-"(12), and dfT = dTjdj etc. 
Equation (|4.20|) becomes 



UJ 



L 2 (h 



T - 



1 — u)\ — nS^u>2 



fd f F-u 2 [V 1 f¥[fd f F + f 2 djf] = 0, (4.24) 



which is the starting point for the evaluation of the two quantities of interest here: the 
autocorrelation exponent A ( |1.3| ) and the late-time scaling form for T . 

For times t\ ^> t 2 both T and / are small. In this limit ( 4.24|) becomes a linear equation 
for T and, following the treatment in ^,11], with the definition Q4.18 ), A can be determined 



as: 



A = d 



H 1 — 10\ — uSqIOi UJq 

^(2) 
^0 



(4.25) 



4fi ,s'. !2i 2 

If one knows u>q, Ux, uj 2 and fi one can determine A. These quantities can be found from an 
analysis of the equal-time correlations, to which we now turn. 

To examine the equal-time order-parameter correlations in the late-time scaling regime 
we set tx = t 2 = t and write ( f4.24| ) in terms of the scaled distance x. To leading order in 
1/Lwe have 



12 



2t r n 1 ~ Ul ~ nS )uJ 2 ... , r . . ^ .... , r ..-2.fi. 



2/2 



(2) 



-fd f T + ^[V./H/c^ + f l d)T\ = 0. (4.26) 



The calculation of the scaling form for T reduces to the solution of the non-linear eigenvalue 
problem (|4.26|) with the eigenvalue /x. The eigenvalue is selected by finding numerically the 
solution of ( |4.26| ) which satisfies the analytically determined boundary behaviour at both 
large and small x. The new aspect to the problem is the presence of the unknowns uq, 
and u)2, a consequence of the incorporation of fluctuations into the model. For n = 2 these 
constants play the role of counter-terms that cancel out the small-x non-analyticities in the 
normalized auxiliary field correlation function /. This procedure fixes uo, iO\ and u>2 in terms 
of fi and d. 

For large x both T and / are small and (|4.26 ) can be linearized. In this regime the 
solution to (|4.26|) is 



T ~ x d - 2X e~ x2 l 2 . 



(4.27) 



The result for the exponent d — 2A appears to be robust. Until now, we have derived results 
valid for arbitrary n > 1. However, the primary goal of this paper is to examine the 0(2) 
model, where there are known qualitative discrepancies with simulation data. With this 
in mind, we now examine the small-x behaviour of the scaling equation (4.26) for the case 
n = 2. For small-x (OB) admits the following general expansion for /: 



f = l + f2X 2 



mx 



\nx 



+ hx A 



1 + ^(1 o 

mx 



lnx 



0(x 6 ). (4.28) 



Non-analyticities appear as a result of the non-zero K2 and K4 coefficients multiplying factors 
of 1/lnx. The non-zero K2 coefficient is particularly important since it is responsible for 
the divergence of the defect-defect correlation function at small x. 

The coefficients of the expansion ( |4.28|) can be determined by examining ( 4.26 ) order- 
by-order at small x. Balancing terms at 0(lnx) gives 



h 
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4/id 



(4.29) 
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This relation between fa and /i is the same one that was found in the original theory [[|. 
This equivalence is a consequence of the simplifications mentioned previously ( |4.19| ) that 
occur for n = 2. At 0(1) we have a equation relating ujq, uj 2 , and K 2 

uo = 2f 2 (l + dK 2 + uj 2 ). (4.30) 

As discussed above, the constant iO\ = — 2 oj 2 for n = 2. If we work with ujq = uj\ = oj 2 = 
then (^4.30|) implies K 2 = —1/d. This is simply the gaussian model examined previously [|J, 
whose non-zero value for K 2 results in the divergent small- x behaviour of the defect-defect 
correlation function g(x). 

Now, however, we can insist that, at 0(x 2 ), f is analytic and enforce K 2 = 0. This 
choice produces the following relation between loq and uj 2 : 

u = 2f 2 (l + uj 2 ). (4.31) 

The small-x divergence in the defect-defect correlation function is now eliminated. The 
leading correction to the 0(x 2 ) term in / is then the 0{x A ) term with a coefficient 

_ d + 3 + u 2 2 f 2 (2 + u ) 
h ~ 2(d + 2) {h) ~ 4(rf + 2) (432) 

We can go further and insist that / is analytic at small-x up to 0(x 4 ). Enforcing K\ = 
allows one to arrive at a complicated expression for uj 2 in terms of \x and d only. 

In section VI we will consider enforcing K 2 = through various choices for the u>e 
and we will numerically solve the associated eigenvalue problem. Before doing this though 
we complete our discussion of the theory by examining in detail the correlations of the 
fluctuations and their relationship to the order-parameter correlations. In the process, we 
will establish the important constraint on the u>£ parameters discussed earlier. 

V. ANALYSIS OF FLUCTUATION CORRELATIONS 



In addition to order parameter correlations the theory also completely describes corre- 
lations in the fluctuation field u. There are two types of equal-time fluctuation correlations 
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that are of interest to us. The first describes cross-correlations between the a and u fields 
and is defined as 

C u0 (l2) = (u(R 1 ,t)-a(R 2 ,t)). (5.1) 

The second describes correlations of the fluctuation field with itself and is given by 

6^(12) = (ui(R 1} t) Uj (R 2 , t)). (5.2) 

As we will see later these quantities are closely related in the scaling regime. One can deduce 
equations of motion for both C u0 and C uu by using the equations of motion ( |3.3|) and ( |3.5|) 
for a and u. For equal-times one has 

^C u0 (12) = V?C u0 (12) - (W y (l)« i (l)o- < (2)) -C eo (12) 

+V^C M0 (12) - C u2 (12) + C ue (12) (5.3) 

and 

-i o 1 1 

--CU12) = V?C UU (12) - -(^(1)^(1)^(2)} - -C eu (12) (5.4) 
2 at n n 

where in the last equation we have used the translation invariance in space to combine two 
equivalent terms. We have also defined 

C eo (12) = (0(1) • a{2)) (5.5) 
C U2 {12) = (3(1) • V 2 m a{2)) (5.6) 
C ue (12) = «l)-e(2)). (5.7) 

We can solve (|5.3| ) and (|5.4|) to determine the correlations for the fluctuation field u if we 
make the additional assumption that u is also a gaussian field. In particular we assume that 
m and u are coupled gaussian fields satisfying 

5 5 
( mi (l)g(m,u)) = C (13)(^^g(m,u)} + C mu (l3)(j^G(m,u)} (5.8) 

and 
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5 5 
( Ul (l)g(m,u)) = C um (13)(— —g(m,u)) + C uu (l3)(——G(m,u)) (5.9) 

dmi{6) 0Ui{6) 

where Q is a general function of fa and u, C um is defined as 

5 l3 C um (l2) = (^(1^(2)), (5.10) 

and integrations over R 3 and £3 are implied. Using these identities all correlation functions 
depending on m and u are determined in terms of Co, C um , and C uu . Thus, we will see that 
equations fl5.3|) and (|5.4|) can be expressed in terms of C u0 , C uu and averages over functions 



of m alone. This is the first step in determining ( |5.1|) and (|5.2|) . The second step is to 
evaluate the averages over m, which can all be expressed in terms of Co, a quantity known 
from our analysis in the last section. The final step is to analyze the equations resulting 
from (|5.3|) and ( |5.4|) in the late-time scaling regime and extract the scaling functions. 

We begin by expressing all of the correlation functions involving a power of u which 
appear in ( |5.3|) and (|5~4] ) in terms of C n o(12), C MU (12) and averages over rh. Using the 
identity ( |5.9| ) we have 

C u0 (12) = C um (12)M 1 (5.11) 
C u2 (12) = C um (12)M 3 (5.12) 
C ue (12) = C um (12)fii (5.13) 
<Wk i (l)u i (l)<T i (2)> = C um (ll)lU 1 (12) + C um (12)W 2 (12) (5.14) 
(^•(1)^(1)^(2)) = C uu (l2)n 2 + C um (21)C um (ll)fi 3 (5.15) 

where we define 

Mi = (V m -*(l)) (5.16) 

M 3 = <V m • V 2 m a{\)) (5.17) 

= -f Mi + £ uj e (V m ■ V^a(l)[V 1 m(l)] 2(£ - 1) ) (5.18) 

n 2 = (Wii(l)) (5.19) 

^3 = (V* m V^(l)> (5.20) 
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Wk(12) = ([V^(l)]a,,(2)> (5.21) 
W 2 (12) = (W ii (l)V^(2)). (5.22) 



We see from ( p. lip that C um can be eliminated in favor of C u0 in equations ( p,12| - |5.l5| ). In 



terms of these various auxiliary functions, equations (|5.ij|) and (|5.4j) become 



^C u0 (12) = V?C M o(12) - ^(HjW^l^Mf 1 - a (12)^ 2 (12)Mr 1 - C eo (12) 

+V^ u0 (12) - a (12)M 3 Mf 1 + a (12)fiiMf 1 (5.23) 



and 



1 s~) T 1 

--CU12) = V?C UM (12) - -C UU (12)Q 2 - -C u0 (12)(C u0 (ll)n 3 Mr 2 + fiiMf 1 ). (5.24) 
lot n n 



The next step is to compute the averages over m ( |5.16| - |5T2~2| ) in the scaling regime. We 



begin with the one-point averages. Except for Q\, which involves spatial gradients, these 
can all be evaluated using 

{A[l7i]) = I j2^^ e ~ X2/(2S0)A[x] - (5 - 25) 
It is straightforward to show that 



and 



*=W*wT<ir (5 - 26) 



Turning to fi 2 we see rather trivially from the form of the stability matrix given by ( f2.9|) 
that 

Q 2 = ql + 0{L- 2 ) (5.28) 

with a InL multiplying the correction term for n = 2. This quantity serves as a mass for the 
longitudinal fluctuations and dominates their determination. For Q% a brief manipulation 
produces the simple result 
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(n-l) 



Q 2 o 



S (t) 



(5.29) 



The last local quantity fix, involves averages over spatial derivatives. The key point in 
handling such quantities is that (mfc(l)Vjmj(l)) = 0. One can then rather easily derive the 
recursion relation, valid for £ > 0, 



((Vm) M G) = nS ( 



(2) 



1 + 



2(£-l) 



nd 



(5.30) 



for a general local function G[m]. With this relation we can evaluate f2i, up to the u>2 term: 

c( 2 ) 



fix = Mi 



UJ 



(5.31) 



L\t) S (t) S (t) 

Turning to the two-point quantities W\ and W2 it is easy to show, using the symmetry 
properties of the order-parameter, that in the scaling regime these reduce to 



and 



Wi(12) = (n-l)ggVo< 



1 



&(1).&(2) 



m(l) 



(5.32) 



^ 2 (12) = g ^ (^[l-(^(l)^(2)) 2 ]). 
m(2) 



(5.33) 



These are new averages to be evaluated. Hereafter, we shall work exclusively with n — 2. 
In this case Wi(12) and W 2 (12) can be evaluated as 



Wx(12) = gg^ 



/ 7T 1 

2^7' 



(5.34) 



while 



W 2 {l2) = ql^. 



(5.35) 



2S 1 + v/W 7, 

We are now in a position to evaluate C u q and C uu in the scaling regime, and to relate 
them to T through Ceo- From the definition (|4.1|) of we see that the scaling ansatz for 
Ceo should have the form 



Ceo = j±F e (x) 



(5.36) 



where, after explicit evaluation, 



TV 



2^ 



(2) K + 2Si 2) LU 2 )fd f F + u 2 [V x ff[fd^ + fd)T\ 



(5.37) 



Looking at the determining equations ( |5.23| ) and ( |5.24| ) we easily see that, as a consequence 
of Q5.36| ), we must take u ~ L~ 2 to leading order. We therefore write 



and 



C u0 (12) = ^F u (x) 



C uu (12) = ^F uu (x). 



With these ansatze (|5.23j ) can be written as 



to leading order, while ( f5.24| ) becomes 



F uu (x) + -F u (0)F u (x) = -\ 
7T qo 
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UJQ 



-{uj 1 + 2S! ) 2) uj 2 ) 



2 F e {x) 



F u {x). 



(5.38) 



(5.39) 



(5.40) 



(5.41) 



2/4 2) 

We see that the quantity F u (0) enters into these equations. If this quantity is to be finite 
then we see that F®(x) can not blow up as x — > 0. Since fdfT does blow up as x — > we 
must choose 



W! + 2Sf ) u 2 = 



(5.42) 



which fixes u\ in terms of ui 2 and tells us, using (pJD , that S Z > = 1 for n = 2, even in the 
presence of these perturbations. We then find that 

2u 2 f 2 — ^0 



^(0) 



2f 2 (l+dK 2 ) 



(5.43) 



where (|4.30| ) has been used. We then have the final results 
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Fu(x) 9 



and 



u [(l + sfl - P)T - /] + 2u 2 f 2 f + u 2 [l + ^l-P][V x f{x)f[fd f F + fd)T\ 

(5.44) 



+ ^F u (0) 

7T 



FJx). 



(5.45) 



Inspection of equations (|5.43|) , ( |5 . 44|) and (|5.45|) show that in the original theory [§J F u (x) = 
F uu (x) = 0, as expected. From the definitions ( |5.2|) and (|5.39|) we must have F uu (0) > 0. In 
the theory with only u Q ^ we have 



Fuu(0) = ^ 1-- 

which is positive. However, if ujq = ( |5.45| ) implies 

2 



(5.46) 



^(0) 



7T 



■[^(0)]' 



(5.47) 



which is negative. Thus within the £ m ax = 2 approximation it is necessary to have ujq ^ 
in order to have a physical theory. For more general G, one must look to the numerical 
solution of ( |4.26| ) to answer the question of the sign of F uu (0). 

Equations (|5.44j ) and ( |5.45| ) explicitly show how correlations in the u field are slaved 
to those of the order-parameter. The universality in ( |5.44| ) and ( |5.45|) is evident, up to 
the non-universal overall factor of l/(?o> which characterizes the flatness of the equilibrium 
minimum in the potential and sets the scale of the fluctuations. 



VI. NUMERICAL ANALYSIS OF THE NON-LINEAR EIGENVALUE PROBLEM 



The eigenvalue problem posed by ( ^4.26| ), subject to the boundary conditions at small- and 
large-a; outlined above, has to be solved numerically. A fourth-order Runge-Kutta integrator 
is used to integrate ( [4.26j ) with initial conditions given by an analytic small-x expansion to 
x = 0.0001. The eigenvalue /i is adjusted until the solution matches onto the gaussian decay 
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( [4.271 ) a t the largest distances. This is now a standard procedure and is essentially the same 
as that used in ||J]. We examine the 0(2) model in two and three spatial dimensions. 

In the original theory || = and the selection of the eigenvalue depended on only 
the conditions outlined above. In that theory K 2 = —l/d, which lead to an unsatisfactory 
non-analyticity in the small-re behaviour of f(x) and ultimately to an unphysical divergence 
in g(x) at small x. In the present theory we can choose 6 so as to eliminate the leading 
non-analyticities and remove the unphysical divergence. Each choice represents a separate 
eigenvalue problem. The simplest choice is to keep only the first term in in by setting 
ujq = 2/2 and ui = for i > 0. For d = 2 the solution to this eigenvalue problem has 
ujq = —1.4620..., while for d = 3 one has u = —1.3450.... Another choice for is to eliminate 
the first term, but keep the next two by setting u = 0, u 1 = 2, u 2 = — 1 and u>£ = for £ > 
2. As mentioned in the previous section, this choice has the unfortunate consequence of 
rendering F uu (0) negative. Both these theories have K 2 = and K4 7^ 0. If we require f(x) 
to be analytic up to 0(x 4 ) we must choose all of u , ui\ and uj 2 to be non-zero, following 
the prescription outlined in Section IV to ensure that K 2 = K4 = 0. The solution to this 
eigenvalue problem has u = —2.6004..., co> 2 = 0.41914... for d = 2 and u = —2.3438..., 
uj 2 = 0.45675... for d = 3. 

Table | contains the eigenvalues /x obtained from these theories. The autocorrelation 
exponents are shown in Table 0. The scaling forms for the order-parameter correlation 
functions T from the various theories are compared in Figure [I] for the 0(2) model in two 
dimensions. The three dimensional results for T are similar, and are not shown. The defect- 
defect correlation function g(x) ( |4. 1 0| ) is calculated from our numerically determined form 
for f(x). The results for the correlations between vortices (n = d = 2) obtained from the 
various theories are compared in Figure 0. Simulation results 0, which are scaled to give 
the best fit to the original theory Q at large x, are shown for comparison. Finally, the 
scaling function F uu ( |5.45| ) is computed using our knowledge of T and /. The results for the 
theory with u ^ and uj 2 7^ are shown in Figure [5] for two and three spatial dimensions. 
For the theory with just ujq 7^ the behaviour of F uu is similar to that shown. 
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VII. DISCUSSION 



It is noteworthy that the scaling solution for we find here for T retains the same qual- 
itative features seen in ||. The differences between the .F's can be attributed mainly to a 
rescaling of the length scale. Since there is always some arbitrariness in the choice of length 
scale when comparing simulation data to theory, we expect that our new results for T will 
fit the simulation data well after rescaling. To our knowledge, simulation results for the 
autocorrelation exponent A of the 0(2) model exist only for two spatial dimensions ||12|| , 
where it is found that A = 1.171. We see from Table [TJ that the original theory || is already 
in excellent agreement with the simulations on this point. It is not surprising then that the 
modified theory makes worse predictions for A than the original theory. The introduction 
of G was not expected to be in any sense a small perturbation. A trend which is counter to 
our expectations is that the discrepancy between the value for A from simulations and the 
value obtained in the theory seems to increase as one includes more terms in 6. It may still 
turn out that the inclusion of higher order terms in I in B does lead to improvement. This 
appears to be a straightforward but tedious calculation. It is also interesting to note that, 
for the theory with ujq = and 002 7^ 0, the prediction for A violates a proposed |20] lower 
bound A > d/2. Despite these problems with the quantitative values obtained for A, there 
is qualitative improvement in g(x) since the small-x divergence for n = d = 2 seen in the 
original theory is removed. The value of g(0) in all the modified theories is too low when 
compared with simulations and this point suggests that some fine-tuning of the theory is 
necessary. Finally, our results for F uu obey the necessary condition F uu (0) > for two of 
our choices of O. We also see that the strength of the fluctuations, characterized by F uu (0), 
increases in lower dimensions, as one might expect. 

We see that the inclusion of fluctuations allows us to render the correlation function / 
of the auxiliary field analytic and we have constructed such a solution up to 0(x 4 ). This, in 
turn, cures the short- distance divergence in g. The theory, at the gaussian level, appears to 
be in better qualitative shape given our development here, however it is at the expense of 
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proper quantitative agreement for the non-equilibrium exponent A. There is evidence f!8 



that the inclusion of post-gaussian corrections lowers the value of A. Thus we hope that 
the tendency for the fluctuations to increase A will be balanced by the introduction of post- 
gaussian terms, resulting in a value for A in reasonable agreement with simulations. We also 
hope that post-gaussian corrections will reduce the magnitude of g(0). It appears that the 
procedure we introduce here leads to a qualitatively more consistent theory. However it is 
also clear that it is unlikely that one can have such a theory and quantitative estimates for 
exponents within the gaussian approximation. One should proceed to look at post-gaussian 
theories. 
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TABLES 



TABLE I. Values for the eij 


^envalue fi from the various theories (in all theories uj\ = 


= -2w 2 ). 


UJ = u 2 = 


ujq / 0, uj 2 = 


CJO = 0, L0 2 = -1 




d = 2 0.94858 


0.53721 


0.76033 


0.42863 


d = 3 0.56837 


0.38931 


0.51434 


n 32544 


TABLE II. Values for the autocorrelation exponent A from the various theories (in 


all theories 


= -2 uj 2 ). 








u = u 2 = 


UJq / 0, w 2 = 


uj = 0, w 2 = -1 


Wo, w 2 / 


d = 2 1.1720 


1.2690 


0.96703 


1.4678 


d = 3 1.6182 


1.6550 


1.4730 


1.7585 
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FIGURES 




FIG. 1. Scaling function J~(x) for the order-parameter correlations in two dimensions. From 
bottom to top, at x = 1, the curves correspond to: the theory with ujq ^ and u>2 / 0; the theory 
with only luq ^ 0; the unmodified theory [Q]; the theory with luq = and UJ2 ^ 0. In all theories 

UJ\ = —2 UJ2- 
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FIG. 2. Scaling function g(x) for the defect-defect correlations in two dimensions. From bottom 
to top at x = the solid curves correspond to: the theory with ujq 7^ and L02 7^ 0; the theory 
with only ujq 7^ 0; the theory with ujq = and u>2 7^ 0) the unmodified theory (diverging). In all 
theories u>± = —2 u>2- The dots represent the simulation results Q for the two-dimensional 0{2) 
model. 
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FIG. 3. Scaling function for the fluctuation correlations for the theory with loq / 0, ll>2 / 
and u>i = —2 u>2. At x = the lower curve is the result for three dimensions and the upper curve 
is the result for two dimensions. 
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